Statistical properties of a localization-delocalization transition induced by correlated 
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The exact probability distributions of the resistance, the conductance and the transmission are 
calculated for the one-dimensional Anderson model with long-range correlated off-diagonal disorder 
at E = 0. It is proved that despite of the Anderson transition in 3D, the functional form of the 
resistance (and its related variables) distribution function does not change when there exists a Metal- 
Insulator transition induced by correlation between disorders. Furthermore, we derive analytically 
all statistical moments of the resistance, the transmission and the Lyapunov Exponent. The growth 
rate of the average and typical resistance decreases when the Hurst exponent H tends to its critical 
value {Her = 1/2) from the insulating regime. In the metallic regime H > 1/2, the distributions 
become independent of size. Therefore, the resistance and the transmission fluctuations do not 
diverge with system size in the thermodynamic limit. 
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I. INTRODUCTION 

According to the pioneering work of Anderson pj, 
any amount of disorder localizes all electrons in one- 
dimensional disordered systems. However, in the case 
of off-diagonal disorder, due to chiral symmetry, there 
is peculiar properties such as, divergence of the density 
of states and the localization length at the band center 
|2jy|jLJ. Scaling studies of the Anderson model with 
diagonal disorder and a number of other models of dis- 
ordered ID systems, have found that the dimensionless 
resistance p = R/T |(| satisfies ln(< p >) oc N and 
< ln(l + p]>oc TV as N the length of the sample goes to 
infinity f^fsBoj . Purely off-diagonal disorder at E = 0, the 
same as diagonal case for < p >, shows an exponential 
growth with the length, but in co ntradic tion with diago- 
nal disorder < ln(l + p) >oc y/N |lflllj |. In the case of 
diagonal disorder, higher moments of the resistance and 
the transmission can be calculated analytically by using 
the generalized transfer matrices method 

In higher dimensions (d > 2), weak disorder does not 
destroy the metallic regime. Only when the strength of 
disorder exceeds a critical value, the electrons become 
localized This phenomena which is called Anderson 
transition, does not depend on the microscopic details of 
the system and is universal. According to the scaling the- 
ory of Abrahams, et al [l^, there exists a single param- 
eter, conductance g, which determines scaling properties 
of g(N). Soon it became clear that the conductance g 
is not a self-averaged quantity. The knowledge of the 
mean value < g > is therefore not sufficient for complete 
description of the transport properties. One has to deal 
with the conductance distribution G(g) 9 14] or equiv- 
alently, with all cumulants of the conductance. This is 
easier in the metallic regime, where G{g) is Gaussian and 
the conductance fluctuations are universal and inde- 
pendent on the value of the mean conductance and/or the 



system size. The width of the distribution depends only 
on the dimension, the physical symmetry of the system 
and depend on the boundary conditions. In the insula- 
tor, large fluctuations of the conductance is characterized 
by the lo g-normal distribution of g. Numerical studies 
proved the system-size invariance of G{g) at 
the critical point, which is consistent with the scaling 
theory of localization. The shape of the distribution is, 
however, not completely understood. 

In addition to the Anderson transition induced by dis- 
order strength, spatially correlation of disorders also can 
change localization behavior of the system. A short-range 
spatially correlation of the onsite disorder causes a dis- 
crete number of extended states in the random dimer 
models. This demonstration revived interests in ID 
disordered models. Special attention has been recently 
paid to the presence of a continuum of extended states 
in ID around th e band c enter in the long-range correlated 
disorder model |5ll9l20| . 

Motivated by the distribution evolution of the Ander- 
son transition, in the present paper, we address to an- 
swer a question about the conductance (and its related 
variables) distribution of the long-range correlated off- 
diagonal disorder along the metallic to insulating transi- 
tion induced by the correlation between disorders. Our 
calculations are done analytically at the band center 
(E = 0). It is proved that despite of the Anderson 
transition in 3D induced by the disorder strength, the 
functional form of the conductance distribution does not 
change when a phase transition occurs by means of dis- 
order correlation. It will be also proved that the conduc- 
tance distribution at the phase transition point and also 
in the metallic regime, has size invariance, while in 3D 
Anderson transition, only the critical distribution is size 
independence. 

Furthermore, we derive analytically all statistical mo- 
ments of the transmission, the resistance and the Lya- 
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punov Exponent and discuss about their fluctuations. It 
has been proved that for long-range correlated hopping 
disorder in the anomalously localized regime H < 1/2, 
< p > and also the typical resistance p — exp[< ln(l + 

\ 1 • • /v 1_2ff /v 1 / 2 -^ 

p) >J — 1 behaves with system size as e and e 

respectively when N — + oo. Here, H refers to the Hurst 
exponent. The critical Hurst exponent where occurs a 
localization-delocalization transition, was proved to be 
H cr — 1/2 in this model [5|- The growth rate of the 
resistance with the length decreases when the Hurst ex- 
ponent goes to the transition point in H cr = 1/2. In 
the extended regime H > 1/2, fluctuations of the trans- 
mission and the resistance do not diverge with system 
size. 

This article is organized as the following sections: Sec- 
tion II introduces our model and our definition of the 
Lyapunov Exponent (L.E.). The transmission distribu- 
tion and its higher moments will be calculated in Section 
III. The conductance distribution and the divergence of 
its moments are discussed in Section IV. Section V fo- 
cuses on the resistance distribution and all its moments. 
Numerical results on the distribution of the L.E. and an- 
alytical results of higher moments at the band center are 
presented in Section VI. Finally, discussions and conclu- 
sions are presented in Section VII. 



energy, the total transfer matrix can be easily derived 
analytically. The function F can be derived as; F = 

Eti lH^) r H%)] where MM =< Hu) >c. a , 

Randomness is imposed on the hit's. By taking a Gaus- 
sian distribution for ln(t)'s, the above sum has a Gaussian 
distribution function with zero mean for enough large 
system size. Therefore, the L.E of such system has a 
semi-Gaussian distribution whose mean is given by: 

V 7T N 

where op can be derived as a function of the pair corre- 
lation function of ln(i)'s[5j. 



2 . ~ 

4 = ~{Ng(0) +2J2(N- i)(-lYg(l)} 

£=1 

0(i-j)=<M~)k(^)>c.«. (4) 
to to 

This equation is converted to the uncorrelated case with 
(g(£) = for I ^ 0) and (g(0) = (rf n{t) ). In this case, the 
variance of F function would be ap — Naf, t y 



II. PURELY HOPPING DISORDER MODEL 



III. TRANSMISSION DISTRIBUTION AT E = 



We consider electrons in ID disordered system within a 
tight binding approximation. The Schroedinger equation 
by the nearest neighbor assumption becomes 



(i) 



where E is the energy corresponding to the electron wave 
function. iV'il 2 is the probability of finding the electron 
at site i, Si are the site potentials which are consid- 
ered here to be zero, and U-u — U^-i = ti the hopping 
terms. Using the transfer matrix method and the above 
Schroedinger equation, the electron wave functions at the 
two ends of the system is related to each other by use of 
total transfer matrix. The total transfer matrix is de- 
fined as:(T/v = rii=i^M-i)' where Tj^_i relates two 
wave functions at neighboring sites. A periodic bound- 
ary condition on hopping terms as t± = tjv+ij causes to 
have unity determinant of total transfer matrix for any 
configuration of disorder. As proved in Ref.p|, the L.E 
with the presence of the above condition can be written 
in the following form. 
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N- 



N 



< IFI >r 



(2) 



where F = ln(|a|) and a is an arbitrary eigenvalue of 
the total transfer matrix. < ... > refers to the average 
of the configurations. In the whole of this paper, our 
calculations are in the special energy, E = 0. In this 



Transmission of a particle through a one-dimensional 
random potential has become a much studied problem 
in the theory of disordered systems. It was shown that 
only the logarithm of the transmission coefficient obeys 
central limit theorem, whereas averages of T and T^ 1 
become unrepresentative of the ensemble for macroscop- 
ically large systems j2^|. The logarithm of the trans- 
mission is proportional to the L.E. and its average scales 
linearly with the length of the system. Transmission fluc- 
tuations are not damped when the sample length goes 
to infinity. On the other hand, its higher moments are 
the same order of its mean. Finite-size L.E. which is 
a self-averaging quantity in the thermodynamic limit, is 
related to the transmission and conductance as the fol- 
lowing form. 



7W = ^KO = 2*Mi + J) (5) 

where g{= T/R) and T are conductance and transmission 
coefficients through the system. Since L.E. distribution 
at the band center has a semi-Gaussian distribution, so 
by use of the above equation, transmission distribution 
t(T) has a log-normal distribution. 



r(T) = 



exp[ 



(InT) 2 



y/2ira F 



T 



(6) 
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The n'th moment of transmission can be easily extracted 
by using the above distribution. 



Higher moments of the transmission in the anoma- 
lously localized region, H < i can be derived as N — > oo. 



< T n >= exp(2n 2 a 2 F )eiic[nJ2a 2 ' 



(7) 



where erfc(a;) is the complementary error function, com- 
monly denoted, is an entire function defined by: 



ei£c(x) = — = 

\/7T 



dt 



(8) 



where for large x, it goes asymptotically to zero as 



S/TTX > 
(1 26 
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and for small x it tends to unity as 
+ ...)). 

Depending on the disorder correlation, the variance of 
(ap) F can be dependent or independent of the length. 
In the present paper, we focus on the long-range correla- 
tion between disorders which results in a phase transition 
from the insulating to metallic phase. For long-range 
correlated disorder, we consider the fluctuations of the 
ln(t)'s are given by the following fl9|: 



.u. 

"t 



<[ln(f)-ln(^)] 2 > c . a .= 2a 1 2 n(t) 



\ l J \2H 



(9) 



where of n ^ is kept fixed for all system sizes p^l2^ | and 
H is Hurst Exponent which determines the strength of 
correlation. The correlation length (£c) is considered 
to be equal to the system size. i,j are the positions of 
the bonds along the chain. The pair correlation function 
arising from Eq.@, results in the following expressions 
for variance of F . 



lim o» F (E = 0)ocl 2 

v I °l»(t) 



1 2H H<h 2 



(10) 



As shown in Ref. 5] , an off-diagonal disordered chain with 
correlation exponent of H < ~ is anomalously localized, 
while in the case of H > |, the system has a metallic 
behavior. Based on the above equation, the transmis- 
sion distribution is independent of size for H > \. Fur- 
thermore, the distribution of all statistical variables in 
the metallic regime of ID off-diagonal disorder (such as 
transmission and conductance) are independent of size, 
while in 3D Anderson model, only the critical distribu- 
tion is size independence. It needs to mention that in 
this model in compared with 3D Anderson model, the 
metallic regime is induced by the long-range correlation 
between disorders. According to Eg . H1U|) . although the 
width of the distributions changes through the phase 
transition, but the whole shape of the distribution func- 
tions does not change. The distributions of the conduc- 
tance and the resistance will be presented in Sections 
(IV, V). 



lim < T" >r 

N >oo 



1 1 

V 2 ™ F n 



<T> 



H < 1/2 (11) 



All moments of the transmission are the same order of its 
mean. According to the Ea. (|10fl for H < |, they go to 

zero with size as a power law behavior (N H ~i). Further- 
more, in the insulating regime, the fractional variance of 
the transmission (the ratio of the variance to the square 
of the mean) in the thermodynamic limit diverges to in- 
finity. 



N 



lim Var{T)/ <T > 2 —> oo 



(12) 



This behavior of the fractional variance shows that < 
T > is not a good representative of the statistical en- 
semble. Therefore the fluctuation of the transmission 
diverges for large samples. 

A perfect transmission through the chain for large 
chain lengths behaves as r(T = 1) ~ N H ~? . The scaling 
behavior of this peak is enough to explain the same be- 
havior for the mean value of the transmission given that 
other peaks of the distribution decay in a similar or faster 
way with size of the chain. In the case of uncorrelated 
disorder which corresponds to H — 0, the mean value of 
the transmission has an inverse square root law in terms 
of the length |24| . This is an interesting behavior since it 
perfectly coincides with scaling predictions for wide wires 
of odd number of transversal modes |2^| . 

In the metallic regime with H > i, the variance of F 
is independent of system size. Therefore, moments of the 
transmission are independent of two variables; size and 
correlation exponent. They only depend on the disorder 
strength. For a weak disorder with condition owt) <C 
1/n, moments of transmission can be expanded in terms 
of the disorder strength. 



< T n > 



N- 



1 



H > 



1 



(13) 



In the metallic regime, also it can be simply shown that 
the fractional variance goes to a constant in the ther- 
modynamic limit. Therefore, the transmission can be a 
representative variable for describing statistical proper- 
ties of the ID metallic system. 



IV. CONDUCTANCE DISTRIBUTION 

The distribution of the conductance can be obtained 
by using its relation with L.E. in Eq.j5]. 



ln (!+i) 2 l 



G(g) = 



1 eX P[ 8^, 

y/2^p~ 5(1+5) 



(14) 
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This distribution converges to a power law form 
l/(g 2 ap), as it should be for large g, but for small g 
the distribution behaves log-normal form as: G(g) cc 

exp[ — " i \/(gap). In the metallic regime, where H > 

i, variance of F will be independent of the system size. 
So, the conductance distribution function will be invari- 
ant as N — » oo. In the insulating regime, H<\, Eq.© 
shows that ap behaves with size as N 1 / 2 ^ 11 in the ther- 
modynamic limit. Therefore, in this regime, the conduc- 
tance distribution for large g, goes off more faster than 
the metallic one. It tends to zero as TV^ -1 / 2 / g 2 . It means 
that in the insulator regime, large conductance occur- 
rence happen rarely with a low probability. For g — ► 0, 
the conductance distribution goes to zero G(g) — ► 0. 
However, two limits of g — > and N — > oo do not com- 
mute with each other. 

Because the conductance distribution decays very slow 
G(g) oc g~ 2 as g — ► oo, all moments of g diverge. Large 
conductance can cause to diverge the mean of g in the av- 
eraging process, although occurrence probability of large 
conductance is very low. Therefore, for omitting the ef- 
fect of large conductance in the mean process, it is better 
to mean the resistance p = 1/g instead of the conduc- 
tance. 



V. RESISTANCE DISTRIBUTION 

The resistance distribution at the band center can be 
also derived by using the definition in Eq.JSJ). 



R(p) 



^/2/na 



1 



(15) 



This distribution goes to (1 — p)/ \j2itap as it should be 
for small resistance p. So, R{p) is an ascending func- 
tion of p when p — > 0. But for large p, this distribution 
converges to the log-normal form. The n'th moment of 
the resistance at E = can be derived as the following 
summation. 



Since ~ e 2n(n i)o> ^ ^ ne f rac ti ona i variance (n=2) 

of the resistance diverges as N — ► oo for H < 1/2. 
Furthermore, all moments of the resistance distribution 
also diverges. According to the results of Anderson 
et al 0, instead of the resistance whose moments di- 
verge, the appropriate quantity is a self-averaged vari- 
able ln(l + p)(= 2Nj) which its fractional variance goes 
to zero as N — > oo. However, the fractional variance of 
the new variable in the hopping disorder model at E = 
tends to a non-zero constant (it/ 2 — 1). This constant 
can be derived by using a relation between the second 
moment and mean of the L.E. as < r y 2 >= n/2 < 7 > 2 . 
Furthermore, the average ln(l+p) has already calculated 
through the Eas. (|3ll0|l as: 



H > h 



< ln(l +p) >oc { (19) 

The typical resistance is introduced as a candidate vari- 
able for experimental resistance. Using the above equa- 
tion, in the thermodynamic limit, this measurable quan- 
tity behaves as p oc e N with the system size for 
H < 1/2. In the special case of H = which corre- 
sponds to the uncorrelated disorder, in good agreement 
with the results of Refs. , the average and typical 

resistances are as e N and . It is clear that the growth 
rate of the average and typical resistance decreases when 
the Hurst exponent tends to its critical value H cr = 1/2. 

In the metallic regime H > 1/2, the resistance only 
grows with the disorder strength. Also, in this regime, 
the fractional variance of the resistance is independent of 
size and also the Hurst exponent. This expression means 
that fluctuations of the resistance and the transmission 
do not increase with size when N — > 00. Therefore, de- 
spite of the insulating regime, they can be considered as 
representative variables for the statistical properties of 
transport. 



VI. LYAPUNOV EXPONENT DISTRIBUTION 



< p n >= E 



£=0 



(-l) l n\ 2 (n-e) 2 * 
(n-l)M\ 



*~i 'exfc[- «/ '2(n- If 'a 2 F 



(16) 

In the anomalously localized regime (H < 1/2), where 
a 2 F goes to infinity in the thermodynamic limit, the above 
equation can be summarized as the following form. 



<p">- (-1)» + 2^f- 



(-1)'"' ,2(n-l?o\ 



(n-l)\t\ 



(17) 



As an immediate consequence of the above general for- 
mula, the average resistance grows with system size as: 



< p >cc e '»(') 



(18) 



With care of our discussions about the transmission, 
the conductance and the resistance properties, the appro- 
priate variable for describing localization properties is the 
Lyapunov Exponent. It can be shown that the L.E. dis- 
tribution function and all its higher moments converge 
for large system sizes. 

Fig.lQJ shows numerical calculations of the L.E. dis- 
tribution for different Hurst Exponents 0. All distribu- 
tion curves have been softened by the Kernel smoothing 
method [2(j without changing any statistical character- 
istic of distributions. It can be seen that the functional 
shape of the L.E. distribution does not change when a 
phase transition occurs by the correlated disorder. The 
pair correlation function arising from Eq.@, results in 
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FIG. 1: Lyapunov Exponent Distribution Function for dif- 
ferent Hurst Exponents at E = 0. The disorder strength is 
considered to be 0.1. The sample size is 10 5 . 



the following expression for the L.E variance as N 



JfT+TR 



H > 



CO. 



(20) 



H < 



This equation shows that correlated disorder causes the 
variance of the L.E. distribution function to converge 
faster than in the uncorrelated case (H=0). It can be 
seen that in good agreement with Anderson et al. |9( , the 
variance of the L.E. in the uncorrelated disorder scales 
according to the law of large numbers as 1/N. 

With care of the semi-Gaussian distribution function 
for the L.E, it can be simply proved that the n'th moment 
of the L.E is proportional to the n'th power of its mean 
as < (7- < 7 >)" >oc< 7 >™. So, Eq.(2Il can be gen- 
eralized to higher moments of the L.E. as the following 
scaling law. 



< (7- < 7 >)" ><x 



q~ln(t) " 



N 



H > 1/2 
H < 1/2 



(21) 



So, all moments of the L.E. in the metallic regime (H > 
1/2) are independent of the Hurst Exponent. Fig.QJ 
shows that the variance of L.E. does not change for H > 
1/2. Of course, at the phase transition point (H = 1/2), 
the width of the distribution will tend to H > 1/2 values 
in the thermodynamic limit. 



VII. CONCLUSION AND DISCUSSION 

The distribution functions of the resistance, the con- 
ductance and the transmission was analytically derived 
for ID off-diagonal Anderson model with long-range cor- 
related disorder at the band center (E=0). Despite of 
the Metal-Insulator transition induced by the disorder 
strength in 3D, the phase transition induced by the cor- 
related disorder in ID does not cause to change the shape 
of the distribution function of the conductance (and all 
related variables). In the metallic regime, the distribu- 
tion functions does not depend on the system size, while 
in 3D Anderson model only the critical distribution is 
size independence. 

We derived analytically all statistical moments of the 
resistance, the transmission and the Lyapunov Exponent. 
In the anomalously localized regime H < 1/2, by means 
of these moments, it was shown that the fluctuations of 
the transmission, the conductance and the resistance di- 
verge in the thermodynamic limit. Convergence of the 
Lyapunov Exponent will happen not only with the sys- 
tem size but with the correlation exponent. The growth 
rate of the average and typical resistances decrease when 
the correlation exponent closes to its transition point 
from insulating to metallic regime. In the metallic regime 
H > 1/2, the resistance and the transmission fluctuations 
do not diverge with the system size. 
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